Asymmetries in the local curvature of the WMAP data. 
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ABSTRACT 

We use the local curvature to investigate the possible existence of non-Gaussianity/asymmetry 
in the WMAP data. Considering the full sky we find results which are consistent with the Gaus- 
sian assumption. However, strong non-Gaussian features emerge when considering the northern 
and southern galactic hemisphere separately, particularly on scales between 1 and 5 degrees. 
Quite interestingly, the maximum non-Gaussianity is found for hemispheres centered near the 
ecliptic poles, which might suggest the presence of some systematic effect. The direction of the 
asymmetry seems consistent with the findings by Eriksen et al. 2004. 

Subject headings: (cosmology:) cosmic microwave background — cosmology: observations — methods: 
data analysis — methods: statistical 



1. Introduction 

Many recent papers (Komatsu et al. 2003; 
Naselsky et al. 2003a,b; Chiang et al. 2003; Park 
2003; Eriksen et al. 2004a; Naselsky et al. 2003c; 
Coles et al. 2003; Copi et al. 2003; Vielva et al. 
2003; Magueijo & Medeiros 2003; Eriksen et al. 
2004b) have tested the Gaussianity of the cos- 



mic microwave background (CMB) on the WMAP 
data. The most striking result from some of these 
papers is the possible existence of an asymmetry 
in the distribution of the CMB fluctuations on 
different hemispheres. For instance (Eriksen et al. 
2004b) find strong non-Gaussianity, as measured 
by the genus, on the northern galactic hemisphere; 
(Vielva et al. 2003) detect, by using a wavelet 
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analysis, a distinct non-Gaussian feature in the 
southern hemisphere; similarly results by (Copi et 
al. 2003; Eriksen et al. 2004a; Park 2003) seem 
inconsistent with the assumption that CMB fluc- 
tuations are isotropic/ Gaussian on various scales. 
The physical implications of these findings are still 
unclear. In this letter we aim at a further inves- 
tigation on this issue by using a different method; 
namely the local curvature approach advocated by 
(Dorc et al. 2003), which we extend from the flat 
sky approximation to the full sky setting (Cabella 
2004). 

The plan of the paper is as follows; in section (2) 
we briefly describe the local curvature formalism 
on the sphere, section (3) outlines our method and 
section (4) collects the results, which we discuss 
in section (5). 

2. Local curvature on the sphere 

Similarly to (Dore et al. 2003) we classify the 
points of the renormalized map: 

in three types: 

• hills where the eigenvalues of the Hessian 
are both positive. 

• lakes where the eigenvalues of the Hessian 
are both negative. 

• saddles where the eigenvalues are of opposite 
sign. 

We can hence probe non-Gaussianity by eval- 
uating the proportions of hills, lakes and saddles 
above a certain threshold i' in the normalized map 
on a given part of the sky (see (Dore et al. 2003)). 
The Hessian on the sphere is calculated using the 
covariant second derivatives (Schmalzing & Gorski 
1998; Eriksen et al. 2004b). The derivatives are 
calculated in spherical harmonic space. More de- 
tails can be found in (Cabella 2004). 

3. The method 

For the analysis performed in this Letter we 
used the publicly available WMAP data^. We 

available at http://lambda.gsfc.nasa.gov 



used the co-added W+V maps and as a consis- 
tency check also the Q1-I-Q2 co-added map. The 
maps where co-added in the way described by 
(Bennett et al. 2003) and the unphysical monopole 
and dipole were removed using the aim coupling 
kernel. In order to lower the noise and to probe 
different scales, we have smoothed the maps with 
Gaussian beams of FWHM in the range 0° to 15°. 
The hills, lakes and saddles were calculated on 512 
simulated maps (in the case where the results are 
less interesting we have limited the number of sim- 
ulations to 128 to save CPU time: this is clearly 
indicated in the tables) as well as on the WMAP 
data using the following procedure: 

1. The map was multiplied with the Kp2 mask 
and normalized. In order to avoid prob- 
lems with smoothing close to the boundaries 
of the mask, the point source holes were 
refilled, mirroring the area just outside of 
the holes. The galaxy was refilled mirroring 
the part just above (for the northern hemi- 
sphere) or below (southern hemisphere). 

2. A spherical harmonic transform was applied 
and the map was smoothed with an ap- 
propriate Gaussian beam in spherical har- 
monic space. An inverse transform created 
the smoothed map; then the mask was ap- 
plied again and a new normalization was per- 
formed. 

3. We calculated the spherical harmonic coeffi- 
cients of this smoothed and masked map in 
order to perform the derivatives in spherical 
harmonic space using a procedure similar to 
(Schmalzing & Gorski 1998; Eriksen et al. 
2004b). 

4. On the derivative maps we applied an ex- 
tended Kp2 sky cut (in order to avoid prob- 
lems with the derivatives close to the bound- 
aries) and for some tests we extended the cut 
even further using the KpO or a 30° x 2 galac- 
tic cut. On these maps we counted the hills, 
lakes and saddles for different levels v in the 
regions of interest. 

Tests on simulations have shown that the bias 
introduced by the refilling procedure for smooth- 
ing around the mask is minimal when compared to 
maps where the smoothing was performed on the 
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full sky. Note also that in the final step, we use 
an extended mask to avoid the unstable deriva- 
tives close to the boundary. The extended mask 
was foiind using the derivative of the mask itself 
(a similar procedure was adopted in (Schmalzing 
& Gorski 1998) to remove the bias introduced by 
smoothing close to the boundaries of a sky cut). 
We included only the pixels in which the deriva- 
tives were above a certain level. The correct level 
was adjusted in such a way that the number count 
of hills, lakes and saddles were not significantly bi- 
ased with respect to simulations with derivatives 
taken on the whole map (the original Kp2 mask re- 
moves 15% of the sky whereas the extended mask 
removes 25%). 

We have hence implemented a statistic: 

U= — 3CT 

where x^, represents the proportions of hills, lakes 
or saddles above a given threshold u and <Tx{i') = 
a/ {x1 — {xi,)"^). We have excluded the highest 
thresholds as there were very few pixels present 
and the statistics becomes very poor at these lev- 
els (see Dore et al. (2003)). 

4. The results 

We applied the test on the full sky W-l-V data 

using the Kp2 mask and a calibration of 128 sim- 
ulations for each smoothing scale. The results are 
reported in table (1), where the numbers repre- 
sent the percentages of simulated maps with a 
higher ■ We see that considering the full sky, 
the maps are consistent with the Gaussian hy- 
pothesis. However, motivated by the recent results 
where non-Gaussianities have been found only on 
the northern galactic hemisphere (Eriksen et al. 
2004b) and with a different method only on the 
southern galactic hemisphere (Vielva et al. 2003), 
we decided to implement the same test on these 
two hemispheres separately. The results are re- 
ported in table (1) for the co- added W-l-V chan- 
nel (512 simulations were used in the calibration, 
except when the percentages are reported as in- 
tegers, in those cases 128 simulations were used). 
The northern hemisphere appears non-Gaussian at 
the 2 — 3(7 level at scales between 1 and 5 degrees. 
The saddles are also non-Gaussian but only at the 



1(7 level and at the smallest scales. The south- 
ern hemisphere shows results consistent with the 
Gaussian hypothesis, except for the scales between 
2 and 4 degrees. Note that for these scales both 
hemispheres taken separately are non-Gaussian, 
whereas considering the full sphere, the results ap- 
pear Gaussian. The reason is that the proportions 
of hills, lakes and saddles are exceptionally high in 
one hemisphere and exceptionally low in the other. 
For the scales between 8 and 15 degrees at which 
the northern hemisphere is Gaussian, the south- 
ern hemisphere produces results with particularly 
small departures from the Gaussian expected val- 
ues; this may be due to a different angular power 
spectrum on the two hemispheres, yielding a lower 
variance than expected in the south (see (Eriksen 
et al. 2004a)). As a consistency check, we per- 
formed the same analysis on the co-added Q1-I-Q2 
map and found very similar results, thus the non- 
Gaussian features do not seem to depend on fre- 
quency channel. To investigate the possibility of a 
contaminating foreground, we have performed the 
same test on the W-l-V channel, using the KpO 
mask and a 30° x 2 galactic cut. The results are 
summarized in table (2). The KpO mask does not 
lower the effect, however with the 30° x 2 galac- 
tic cut, the significance is reduced to the Icr level. 
This does not imply that the feature which caused 
a significant non-Gaussian detection using a larger 
part of the sky is absent outside of this huge cut. 
We will sec at the end of this section that the posi- 
tional dependence of the does not change even 
when using the 30° x 2 galactic mask. Finally we 
also point out that the hill, lake and saddle counts 
are power spectrum dependent. However, trying 
small variations of the model power spectrum in 
the simulations does not seem to change the re- 
sults significantly. 

In order to investigate the exact location of the 
non-Gaussian features, we implemented the tests 
on hemispheres centered at 164 different positions. 
In table (2) we show the results on the maximally 
asymmetric directions. In each simulated map we 
estimated the on each of the 164 differently ori- 
entated hemispheres and we recorded the highest 
X^ value. We did the same for the WMAP W+V 
map with the Kp2 sky cut. Taking into account 
the rotations, we see that the significance of the 
asymmetry is rciduccid, but still at the 2cr level for 
the 3 degree smoothing scale. 
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In figure (1) we represent the dependence of 
non-Gaussian features from different directions for 
the 1 degree smoothing scale using the Kp2 cut. 
The discs represent the centers of the hemispheres; 
the numbers reported indicate the percentages of 
Gaussian simulations which yielded a higher 
value for the lakes for that given hemisphere. Note 
that there are two non-Gaussian areas, both close 
to the ecliptic poles, which could suggest the pres- 
ence of a systematic effect at this scale. For the 
3 degree smoothing scale (not shown), the map 
appears very similar, the main difference being 
that the two non-Gaussian 'poles' are now slightly 
closer to the galactic poles. However for the 5 de- 
gree smoothing scale, the non-Gaussian area in the 
southern hemisphere vanishes. Using the 30° x 2 
galactic cut, the appearance of the map hardly 
changes, suggesting that the same features are still 
present, even if they are only at the Icr signifi- 
cance level. Furthermore, the maps of the for 
the hills have a very similar appearance to those of 
the lakes; also, using the co-added Q1-I-Q2 chan- 
nels, the same structure emerges. Figure (1) is 
strikingly similar to an analogous one in (Eriksen 
et al. 2004a) which shows the distribution of the 
power spectrum and the three-point correlations 
function on the sky. 

5. Conclusions 

In this letter we have further investigated the 
possible existence of non-Gaussianity/asymmetry 
in the WMAP data. Considering the full sky 
data (using an extended Kp2 galaxy cut and point 
source mask) we find results which are consistent 
with the Gaussian assumption. However, strong 
non-Gaussian features (2-3cr) emerge when con- 
sidering the northern and southern galactic hemi- 
sphere separately, particularly on scales between 
1 and 5 degrees where the proportions of hills and 
lakes are extermely high on one hemisphere and 
extremely low on the other. Quite interestingly, 
the maximum asymmetry is found for hemispheres 
centered close to the ecliptic poles, which could 
suggest the presence of a systematic effect. For 
the 5 to 10 degree FWHM smoothing scale, the 
southern hemisphere yields results so close to the 
Gaussian expected values to suspect ovcrcstima- 
tion of the variance in this part of the sky. The 
position of the non-Gaussianity on the sky seems 
very consistent with the asymmetric distribution 



of the power spectrum and the three-point corre- 
lation functions found by (Eriksen et al. 2004a). 
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Fig. 1. — The discs show the centeres of hemi- 
spheres on which the been calculated for 
the lake counts. The numbers indicate the per- 
centage of simulations with a higher for the 
given hemisphere. In the analysis, the extended 
Kp2 cut was used with a smoothing scale of 1° 
FWHM. The two small dots indicate the position 
of the ecliptic poles. 



This 2-column preprint was prepared with the AAS lAT^rjX 
macros v5.2. 
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Table 1 
Results 



FWHM 


hills 


lakes 


saddles 


full sky (Kp2) 


no sm. 


72 


56 


63 


0.5° 


2 


45 


1 


1° 


39 


39 


40 


3° 


24 


16 


21 


5° 


20 


29 


62 


8° 


48 


41 


41 


10° 


56 


83 


55 


15° 


34 


76 


26 


northern gal. hemisphere 


(Kp2) 


no sm. 


3 


8 


16 


1° 


0.8 


6.3 


10.7 


2° 


1 


4 


5 


3° 


0.8 


1.4 


13.5 


4° 


2 


1 


53 


5° 


1.4 


1.8 


60.9 


8° 


10.4 


18.9 


49.4 


10° 


13 


34 


37 


southern gal. hemisphere 


(Kp2) 


no sm. 


28 


40 


98 


1° 


28.1 


37.5 


75.0 


2° 


2 


20 


63 


3° 


2.5 


3.9 


79.1 


4° 


25 


2 


23 


5° 


70.1 


32.6 


56.3 


8° 


92.0 


34.6 


78.3 


10° 


91 


83 


99 
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Table 2 
Results 



FWHM hills lakes saddles 



northern gal. hemisphere (KpO) 

P To 7^8 16^8 

3° 1.6 1.2 20.9 

5° 3.1 1.2 46.3 

8° 13 15 50 
southern gal. hemisphere (KpO) 

■p 32^0 504 73^2 

3° 2.9 7.8 80.9 

5° 58.4 50.0 93.2 

8° 77 72 91 

northern gal. hemisphere (30° x 2 gal. cut) 

■p 13 30 38 

3° 14.6 9.6 41.0 

5° 14.5 7.4 54.5 

8° 23 45 45 

southern gal. hemisphere (30° x 2 gal. cut) 

■p 49 65 94 

3° 8.6 22.1 89.0 

5° 68.8 74.4 100.0 

8° 68 84 91 
maximally non-Gaussian hemisphere (Kp2) 

P Ol iJl 30 

3° 5.7 4.5 26.2 

5° 10.4 9.4 59.8 
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